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ABSTRACT 



^ The results of [1] are used to obtain full asymptotic expansions of Feynman diagrams 



renormalized within the MS-scheme in the regimes when some of the masses and exter- 
nal momenta are large with respect to the others. The large momenta are Euclidean, 
and the expanded diagrams are regarded as distributions with respect to them. The 
small masses may be equal to zero. The As-operation for integrals is defined and a 
simple combinatorial techniques is developed to study its exponentiation. The As- 
operation is used to obtain the corresponding expansions of arbitrary Green func- 
tions. Such expansions generalize and improve upon the well-known short- distance 
operator-product expansions, the decoupling theorem etc.; e.g. the low-energy effective 
Lagrangians are obtained to all orders of the inverse heavy mass. The obtained expan- 
sions possess the property of perfect factorization of large and small parameters, which 
is essential for meaningful applications to phenomenology. As an auxiliary tool, the 
inversion of the -R-operation is constructed. The results are valid for arbitrary QFT 
models. 



1 on leave from Institute for Nuclear Research of the Russian Academy of Sciences, Moscow 117312, 
Russia 
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1 Introduction 

In the preceding paper [1] the problem of obtaining asymptotic expansions of multiloop 
Feynman diagrams was analyzed. From the point of view of phenomenology, it was 
stressed that such expansions must possess the property of "perfect factorization" of 
large and small dimensional parameters. From the point of view of mathematics, it 
was shown that the key technical problem ( "the Master problem" ) is that of obtaining 
asymptotic expansions — in powers and logarithms of the expansions parameter — of 
products of singular functions in the sense of distributions. An explicit solution for the 
Master problem was obtained in the form of the so-called As-operation for products 
of singular functions, which acts on Euclidean momentum-space Feynman integrands 
and yields an expansion whose all terms are well-defined distributions (containing, in 
particular, 5-functional contributions with coefficients depending non-analytically on 
the expansion parameter) so that the expansion allows termwise integrations with test 
functions. 

The purpose of the present paper is to use the results of [1] and obtain explicit 
expansions of multiloop Feynman integrals in the form of As-operation for integrals. 1 
It turns out that the combinatorial structure of the As-operation for integrals is very 
similar to that of the i?-operation 2 . This similarity allows one to derive expansions of 
perturbative Green functions in a global form with the help of the same techniques as 
used in studying the exponentiation of the i?-operation. 

The present paper is a revised and simplified version of our publications [2] , [3] and 
[4] which have been widely discussed in the literature [20]-[22]. 3 The revision has not 
affected the results but only the order of presentation. The simplification is due to the 
fact that the starting point of [2] was the so-called EA-expansion for UV-renormalized 
integrals, which is a composition of i?-operation and the As-operation in the sense 
of the present paper. Correspondingly, the As-operation as defined in [2] mixed up 
UV-renormalization and the expansion algorithm. In the present version we define the 
As-operation and study its exponentiation directly for the unrenormalized integrals, 
which results in further simplifications. The attractiveness of the As-operation is due 
to the fact that it serves as a combinatorial organizing center of the theory and fully 
exhibits the structure of the expansions, which is invaluable for applications where 
diagram-by-diagram analysis is necessary to enhance reliability of calculations. Its def- 
inition is also remarkably simple within dimensional regularization, which trivializes 
the study of exponentiation. It should be stressed, however, that the use of dimen- 
sional regularization is by no means essential for the expansion problem; a separate 
publication is devoted to a study of this point — cf. [14], [15], and especially [17]. In the 
latter work Euclidean asymptotic expansions were studied in regularization indepen- 

1 The two ^-operations are different representations of essentially the same operation, which 
justifies using the same name to denote them: the As-operation for integrals is an integrated version 
of the As-operation for products — cf. subsect. 5.1 below. 

2 Note that the As-operation for products also has close parallels with the i?-operation — but with 
the i?-operation in position representation. 

3 OPE in the MS-scheme was also studied in [19]. 
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dent manner within the framework of generalized MS-schemes [15]. The main result of 
[17] is an expansion formula whose combinatorial structure is practically identical to 
that of the As-operation of the present paper. Therefore, all our combinatorial results 
concerning exponentiation of the As-operation are immediately applicable to [17]. 

The plan of the present paper is as follows. It consists of five sections, each section 
is subdivided into subsections and contains a prologue where further information can 
be found. In the first three sections various aspects of combinatorial structure of the 
Reoperation in the MS-scheme are studied. In sect. 2 the exponentiation formula for 
the i?-operation is obtained using the functional techniques. The result is, of course, 
well-known. However, since we wish to use similar arguments in studying the As- 
operation, we feel it is useful to present an explicit derivation in such a way as to make 
it immediately applicable to the case of the As-operation studied in sects. 4 and 5. In 
sect. 3 the inversion of the Reoperation is derived, and in sect. 4 the renormalization 
of Green functions with multiple local operator insertions is considered. 

In sects. 4 and 5 we study the combinatorial structure of asymptotic expansions of 
Feynman diagrams. In sect. 5 we consider the special case of heavy mass expansions. 
The As-operation for integrals is introduced and its exponentiation is proved. Using 
the techniques developed in sects. 1-3, the effective low-energy Lagrangian is presented 
in an explicitly convergent form to all orders in the inverse heavy mass. (It should be 
stressed that in this paper we only deal with purely combinatorial aspects of the theory. 
A fuller justification of a technical assumption (see subsect. 5.8 below) is presented in 
[17] (for an informal discussion see [16].) In sect. 6 the results of sect. 5 are extended 
to comprise the case of large external momenta, so that the familiar short- distance 
operator product expansion is reproduced and its generalizations are obtained. 

Note that our final formulae immediately provide the very useful explicit expressions 
of OPE coefficient functions first announced in [11] (explicit examples of calculations 
going as far as 3-loop approximation can be found in [18]; for further references see 
also [23]). 

The notations are on the whole consistent with those of [1]. 
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COMBINATORICS OF THE ^-OPERATION 



IN THE MS-SCHEME 



2 /^-operation in the MS-scheme 



In this section we study the combinatorial structure of the i?-operation applied to the 
perturbation series as a whole. The fact that the i?-operation is equivalent to adding 
counterterms to the Lagrangian is of course well-known. But we are interested in 
exact formulae expressing this fact that would be convenient for practical calculations 
and valid for the MS-scheme of UV renormalization. Therefore we reproduce a simple 
proof of the exponentiation of the -R-operation in the MS-scheme — see (2.25) below — as 
given in [2] . Our second aim is to present the derivation of the exponentiation formula 
(2.25) in a form sufficiently general to make it applicable to the case of the As-operation 
studied in sects. 4 and 5. To make the paper more accessible to practitioners of applied 
QFT, our presentation is rather detailed. 

We start in subsect. 2.1 with a description of the i?-operation in the MS-scheme. 
Our definitions are equivalent to the standard ones [5], and the differences are mainly 
notational. In subsect. 2.2 a convenient order of enumeration of UV-subgraphs is 
fixed in preparation to studying the global structure of the i?-operation. Subsect. 2.3 
summarizes some results of the functional techniques that will be needed (for more 
details see e.g. [7]). In subsect. 2.4 a formal expression of the statement that the R- 
operation is equivalent to adding counterterms to the Lagrangian, the exponentiation 
formula (2.25), is derived; our reasoning is an extension of the arguments of [6] — for 
comments and comparison see subsect. 2.4). In subsect. 2.5 the results are extended 
to Green functions of fields. Complications due to composite operator insertions will 
be discussed in sect. 4. 

2.1 Basic definition of the /^-operation 

The effect of the i?-operation in the MS-scheme on a single Feynman integral G can 
be expressed by the formula 



Let us explain it. 

One associates a Feynman graph, that will be referred to as "the graph G" , with 
the integral G, and vice versa. The correspondence between graphs and integrals is 
established by Feynman rules. 

Let us fix how /i, the 't Hooft unit of mass, appears in the integrals. In momentum 
representation, having taken into account momentum conservation in all the vertices of 




(2.1) 
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the graph, one can obtain the integrand of the corresponding unrenormalized Feynman 
integral. Then for each internal momentum integration one replaces the measure as 



d A p 2e d D p 



■ A* 7— tf, 2e = A-D, (2.2) 
(2vr) 4 P (2tt) d ' ' V ; 

and assumes that integrations (as well as the Dirac algebra etc.) are to be done within 
dimensional regularization. The role of fi is to preserve the dimensionality of the 
integral whole under regularization. 

An UV-subgraph g of G consists of some vertices of the graph G and some of its 
lines attached by both ends to the vertices of g. If g consists of one vertex and no 
lines then we call it elementary vertex. If g has at least one loop and is connected 
and one-particle-irreducible then we call g a proper UV-subgraph. (Note that if G is 
one-particle-irreducible then one of its proper UV-subgraphs coincides with G.) 

An V V -partition of G is a set of its UV-subgraphs {g a } such that each vertex of 
G belongs to one and only one g a . Note that one of the UV-partitions consists of 
elementary vertices only, and another one consists of only one subgraph, namely, G as 
a whole. Summation in (2.1) runs over all UV-partitions of G. 

Now we have to define the operation A\jy°g where g is an UV-subgraph of G. By 
definition, A\jy°g = 9 if 9 is an elementary vertex, and A\jy°g = if g is not an UV- 
subgraph. If g is a proper UV-subgraph, then the action of A uv on g can be described 
as follows. The graph for (Auv°flO x (G/g) can be obtained from the graph for G 
by shrinking g in G to a point. With the new vertex thus obtained one associates a 
factor — the UV-counterterm — that is obtained from the integral corresponding to g via 
a special algorithm which need not be specified here (see, however, below). Using the 
Feynman rules thus extended one finally builds up the integral for (A uv °(?) x (G/g). 
Note that the above rule for \x is still operative, i.e. the number of the factors jj? e is 
equal to the number of loops of G with g shrunk to a point. 

For the sake of completeness let us describe the recipe for evaluating the UV- 
counterterm for a given proper UV-subgraph g. If one rewrites the i?-operation as 

R°<? = R'og + A uv °</, (2.3) 

then Ajjy°g is precisely the UV-counterterm to be evaluated while R' involves only 
counterterms for subgraphs of g which are assumed to be known already. Since we are 
working within dimensional regularization, all our integrals depend parametrically on 
D, the complex-valued dimension of space-time. Let K be the operation that picks out 
the pole part at D = 4 of any function on which it acts. R°g should be finite at D — 4, 
therefore K°R°g = 0. On the other hand, the UV-counterterms in the MS-scheme are 
pure poles. So, applying K to (2.3), we get: 

A uv °# = -K°R'°g. (2.4) 



Eqs. (2.1), (2.3) and (2.4) provide a convenient description of the /^-operation in the 
MS-scheme. We would like to stress that all the manipulations are to be done in 
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momentum representation, and that the rule for // remains operative. Also recall that 
(2.4) is a polynomial of masses and external momenta of g and is independent of //. 

The most important points for us here are: 

(i) A\jy°g, whatever the subgraph g is, depends only on g but not on G as a whole. 

(ii) Summation in (2.1) runs over all partitions of G. 

(Hi) An UV-subgraph is defined as a subset of vertices and some (not necessarily 
all) of the lines connecting the vertices of this subset. 

On the other hand, the fact that Auv is non-trivial only on proper UV-subgraphs 
will be completely irrelevant. 



2.2 Enumerating UV-subgraphs 

Let us transform (2.1) to a more convenient form. Indeed, there are UV-subgraphs 
that have the same set of vertices but differ only with respect of the number and/or 
arrangement of lines connecting them (e.g. when a line of a subgraph can be removed 
without changing its status as 1PI) — but not vice versa: if the sets of lines of two 
subgraphs are identical, than their sets of vertices are identical, too. This allows one 
to perfrom enumeration of UV-subgraphs in (2.1) in two stages: first, one enumerates 
sets of vertices, then, sets of lines connecting those vertices. 

Indeed, let v be a subset of vertices of G and 

A v oG=J2(A VY og) x (G/g), (2.5) 

where summation runs over all UV-subgraphs g with the same set of vertices, v. De- 
noting the set of all vertices of G as V, we can rewrite (2.1) as: 

RoG= J2 (U A v a )°G. (2.6) 

v=Uv a \ a / 

D a ni>/3=0, for aj^fl 

2.3 T-products and functional techniques 

Let 0(x) be a free field operator; (p(x) will denote the corresponding classical field. Let 
F(ip) be a functional of the form 

F (f) = J2 F n (xi . . . x n )ip(xi) . . . ip(x n )dxi . . . dx n . (2.7) 

n J 

Replacing the product of </?'s in (2.7) by the T-product of </?'s one obtains an operator 
which we denote as TF(<p). The Wick theorem allows one to reexpand TF(ip) in terms 
of the Wick normal products of (p : 

TF($) = I F^( Xl . . . x n ) N [<p( Xl ) . . . <p(x n )} dx ± ... dx n , (2.8) 

n J 
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where are related to F n in a certain way (see below). 

The generating functional for the coefficient functions is defined as: 

^V) = E / F n(%l ■ ■ ■ XnMx!) . . . if{x n )dx 1 . . . dx n . (2.9) 

The relation between F(</?) and F N (ip) is expressed, by definition, by 

TF(<p) = N[F N (^)]^. (2.10) 

The well-known formal expression of the Wick theorem due to Khori (see [7] and refs. 
therein) reads: 

F»=e<F(y,), (2.11) 



where 



and 



I = - [ dxdy-^—A c (x : y)-^— = -5A C 5, (2.12) 



A c (x,y)=<T0(x)^y)>o. (2.13) 

(Extension of all our formulae to the most general case of fermionic and complex fields 
is straightforward, therefore our analysis will be quite general.) 

Our purpose is to study Green functions, but first it is convenient to consider the 
(off-shell) S'-matrix. The unrenormalized ^-matrix is expressed as: 

S(y3)=Texp[L(y3)], (2.14) 

where L(ip) denotes the interaction Lagrangian integrated over space-time and multi- 
plied by i (the normal ordering of fields in L(<p) is not assumed). 

The generating functional S N is defined as (cf. (2.10)): 

Texp[L(y3)] = N[S N M]^ (2.15) 

and from (2.11): 

S N = e l exp[L(<p)]. (2.16) 

Again we have to set the rules for /x, in order to establish connection between our earlier 
definitions of Feynman integrals and the expressions generated using the functional 
technique. We assume that: 

(i) all momentum and space-time integrations, 7-matrices etc. are D-dimensional; 

(ii) the dimensions of the fields if and the coupling constants are always the canon- 
ical 4-dimensional ones; 

(Hi) each propagator bears an extra factor /i 2e ; note that (2.12) is dimensionless 
since we assume that 8/8ip(x) x ip(y) = 5(x — y) with D-dimensional 5-function; 

(iv) the interaction Lagrangian L(tp) in (2.14) contains the factor fi~ 2e and is 
therefore dimensionless; note that condition (Hi) is equivalent to (iv) extended to 
the quadratic part of the full Lagrangian ( cf. (2.35) below). 
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Given these rules, a connected Feynman integral will have the factor /i 2 ^' -1 -* where 
/ is the number of its internal momentum integrations (loops) — which differs by fi~ 2e 
from what was postulated in subsect. 2.1. Such extra fi~ 2e to each connected component 
can be easily taken into account in what follows. 



2.4 Exponentiation of the /^-operation 



The renormalized S-matrix is obtained by applying the i?-operation to each Feynman 
integral contributing to the coefficient functions of S N . Thus, the starting point of our 
analysis is the following expression: 



OO 1 

S R M = R°[T exp(L(y.))] = £ — RoTL^), 

N=0 iV • 

where we have introduced the convenient notation 



(2.17) 



(2.18) 



which will be used systematically in the context of the functional techniques. 



Evaluating functional derivatives, one obtains a sum of terms, each one having a 
graphic representation. In the resulting expression one only has to replace products 
of the classical fields y? by the normal (N-) products of the free fields ip, in order to 
obtain the operator of the S-matrix. (Note that since we wish to use the MS-scheme 
for UV renormalization, the S'-matrix will not satisfy correct normalization conditions 
automatically, but this is of no importance because our study of the S'-matrix is only 
an intermediate step in the study of Green functions.) 

Note that in (2.17) each L(cp) corresponds to a vertex while each / (see (2.12)) 
generates one line of a Feynman graph. This allows one to conveniently perform all the 
enumerations of UV-partitions inherent in the definition of the i?-operation as given 
by eqs. (2.1)-(2.4). 

First one has to enumerate all possible partitions of N vertices into groups of vertices 
of various sizes (cf. (2.6)). N vertices can be split into n\ groups of 1 vertex, n 2 groups 
of 2 vertices . . . groups of N vertices in 



' N ' 
n x . . . n N , 



■ N 
k=l 



-1 



J2 k x n k = N, 

k 



(2.19) 



ways. Therefore, applying (2.6) to (2.17) we get: 



OO 1 

tv -L y 

N=0 ' ™i--- n jv 




. . . (a°l n ) 



n N 



(2.20) 



(Note the abuse of notation: A acts on subgraphs consisting of vertices and lines — but 
the latter will only be generated by T. However, we use such a representation only for 
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the purposes of combinatorial enumeration so that no problems should arise.) Now, 
for each group of vertices we have to expand A according to (2.5). We have to separate 
the lines connecting vertices inside each group A°L k from the lines connecting vertices 
of different groups, and among the former to enumerate all subsets of lines that will 
form UV-subgraphs on which Auv acts. This can be done by ascribing different labels 
to the n\ + . . . + tin groups of vertices and marking the field f in each group by this 
label. Then we have: 



exp(-5A c 5) 



1 



exp | ^ A ^/3 + E ^A c 8 a ] I]A^(^) 

a<(3 a 1 



exp ( ]T 5 a A c 5 p n A ° (exp(-5 a A c 5 a L k "( Va ) 



ip a =(p 



<p a =<p 



(2.21) 



(here 5 a = S/5ip a while the index k a takes the value 1 rii times, . . . , the value N, tin 
times). In (2.21), each operator l a = 5 a A c 5 a generates one line attached to vertices of 
the a-th group. There are Q£) = M!(M'!(M - M'y.)- 1 ways to choose M' lines from 
a set of M lines. Therefore, 



Me l F(<p))= E m\HI M F{v)) 



M=0 



M=0 1V1 ■ M'=0 




M-M' 



(2.22) 



where 



^uv 



H " 2e A uv /i 2e . (2.23) 

The origin of the seemingly bizarre /x-factors in (2.23) is as follows. The i?-operation 
and the operator Auv should be applied to "standard" dimensionally regularized Feyn- 
man integrals containing one factor fi 2e per each loop. But the expression e l F(ip) gener- 
ates integrals with one such factor lacking per each connected component — cf. the end 
of subsect. 2.3. To remedy this, the powers of \x are introduced into (2.23), where it has 
been taken into account that Auv is non-trivial only on integrals with one connected 
component. It should be noted that in the context of the functional techniques Auv is 
systematically replaced by A uv , so that the superscript / can and will henceforth be 
omitted. 



Using (2.22), we rewrite (2.21) as: 

eq.(2.21) = TnA UV °(Ti>(v?)), 



(2.24) 



where each T, as before, is the functional differential operator defined by (2.18). Sub- 
stituting this into (2.20), we finally obtain: 



S R (<p) =Texp[L R ], 



(2.25) 



2 R-OPERATION IN THE MS-SCHEME 



10 



where 

L R = A uv »(Te^)-l). (2.26) 

Note that A\jy°(L(ip)) = L(tp) (recall that Auv is a unit operation on elementary 
vertices), therefore Lr can be represented as: 

L R = L((p) + "divergent UV-counterterms" , (2.27) 

where 

"divergent UV-counterterms" = A uv o [Te L(v) - L(tp) - l] (2.28) 

(note that each term in the bracketed expression on the r.h.s. of (2.28) that gives a non- 
zero result after application of A uv has at least one loop). So, indeed, the effect of the 
i?-operation is equivalent to adding some counterterms to the interaction Lagrangian. 
Eq.(2.28) provides a convenient explicit expression for them. 

For clarity's sake, let us explain the algorithm encoded in (2.28): 

(1) one writes down all the contributions to T exp[L(y?)] = exp(|5A c £) exp[L(</?)]; 

(2) then one discards all the terms except those that correspond to connected 1PI 
graphs with at least one loop; 

(3) each of the terms left has the form: 

^ J (II d Pi )(2n) D i5(J2Pi)FI( Pl . ..p k )<p(pi) • ..<p(p k ), (2.29) 

i 

where 

(p{p) = (2ix)- D I e tpx d D x<p(x), (2.30) 

and S(J2Pi) expresses momentum conservation while FI(jpi . . .pk) is a connected 1PI 
Feynman integral constructed according to the rules described after (2.1). To FI, one 
should apply Auv defined after (2.2), and Auv replaces FI in (2.29) by a polynomial 
of pi with coefficients that are divergent in the limit D — > 4. That is, Auv transforms 
(2.29) into an integral over x of local products of (p(x) and its derivatives. (Note that 
in the case of the vacuum graphs (k = 0) the ^-function degenerates into an ill-defined 
factor 6(0), and a special axiom should fix a recipe for handling it. Such a recipe is 
completely non-interfering with what we are doing: the integral FI is correctly defined 
even in this case for both the R- and As-operations to yield meaningful results when 
applied to it.) 

The method for resolving combinatorics of the i?-operation that we have used is 
an extension to the more complicated case of the MS-scheme of a reasoning from 
subsect. IV. 1.3 of [6]. In [6] the momentum subtraction scheme [8] was used (so that 
only one term on the r.h.s. of (2.4) was retained, namely, the one corresponding to 
the UV-subgraph containing all the lines of the graph G which connect the vertices 
from the subset v), and it was assumed that the Lagrangian is normally ordered. As 
a result, in [6] the counterterms in (2.28) were to be understood as normal products, 
while no such normal ordering is assumed in our case. 
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2.5 /^-operation on Green functions 



Let us now turn to Green functions. To obtain the generating functional G(J) for the 
Green functions one should add a source term to the Lagrangian: 



L — > L + ipj, 

and evaluate the vacuum average of the resulting S'-matrix. 



(2.31) 



The term (pj generates a vertex that can be either isolated from, or connected by 
only one line with the rest of the graph. In either case the resulting graph is nullified 
by Auv unless it consists of a single vertex ipj or does not contain such vertices at all. 
Formally: 

A uv ° (T e L+ipJ ) = A uv ° (Te L ) + A uv °(y? J) = A uv ° (Te L ) + ipj. (2.32) 

The last equation is due to the fact that Auv does not affect single vertices. Therefore, 
for the generating functional of the MS-renormalized Green functions one has: 



G(J) = R° < T exp[L(^) + <pj] > = (T exp[L R + pJ])^ =Q . 
It may be helpful to transform (2.33) to the form of a functional integral: 



(2.33) 



1 



(T exp[xM + y?J]L =0 = exp(-<5A c <5) exp 



X {iJ 



= exp 



exp ( -^A C J 



exp 



Mm. 



exp [tpj] J 
dipexp[L {ree + ipj] 



where 



= J d(pexp[L {rcc + x(f) + V-J], 



(2.34) 



(2.35) 



is the free classical action (multiplied by ifi 2e ) of the field ip. Denoting by L tot the full 
classical action (multiplied by ifi~ 2e ) of the model: 

Ltot = L irec + L, (2.36) 

and using (2.27), we have for (2.33): 

G(J) = J dip exp[L to t(p) + "divergent UV-counterterms" + pJ]. (2.37) 

This equation together with (2.28) provides a convenient link between practical cal- 
culations of counterterms and the analysis of Green functions by the renormalization 
group method. 
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3 Inversion of the /^-operation and the ^-mapping 

Since the effect of the i?-operation consists in adding to a graph G a linear combination 
of graphs with a lesser number of loops, it turns out possible to construct an inversion of 
the Reoperation — in very much the same way as it is always possible to invert a triangle 
matrix with units on the main diagonal. We construct the inverted Reoperation, and to 
this end the so-called formalism of ^-mapping is developed. This formalism proves to 
be a useful conceptual tool for studying various aspects of renormalization owing to the 
universality of the ^-mapping, i.e. owing to the fact that it contains full information on 
the Reoperation in a given subtraction scheme for arbitrary models with a given field 
content. The main application of the inverted Reoperation will be the construction 
of Euclidean asymptotic expansion of Feynman diagrams in an explicitly convergent 
form. 

After introducing some definitions in subsect. 3.1, in subsect. 3.2 the ^-mapping is 
defined. In subsects. 3.3-3.5 its inversion, the mapping £ _1 is proved to exist. A useful 
property of is pointed out in subsect. 3.6, and in subsect. 3.7 we briefly discuss 
how the structure of the renormalization group emerges from the point of view of the 
formalism of ^-mappings. In subsect. 3.8 the inverted Reoperation is defined and an 
explicit algorithm for calculating the corresponding counterterms is derived. 



3.1 Local operators, local functionals, and Lagrangian func- 
tionals 

Let jk{x) be a full linearly independent set of local products of (p(x) and its derivatives. 
In momentum representation, one has: 



3k 



(q) = J d D Pl . . . d D p m (2tt) V 2e 5 (q + f>*) ^ 



X<p(pi)...<p(Pm), (3.1) 

with (p{p) defined in (2.30); Vk is a polynomial of p\ . . .p m . For convenience we assume 
that Vk are polynomials of the masses of the model and are allowed to contain factors 
depending only on q = — Vk — or in terms of the coordinate representation, the local 
operators jk(x) are allowed to contain full derivatives in x. 

We say that £ is a local functional of fields if 

£ = £ / dqC n {q)j n {q), (3.2) 

where C n (q) are coefficient functions that are independent of the fields. 
We call £ defined in (3.2) a Lagrangian functional if 

= 9nS{q), (3.3) 

where g n are coupling constants independent of q. For example, L and L R in (2.14)- 
(2.16) and (2.26) are Lagrangian functionals. 
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3.2 The £- mapping 

One can see that the proof of eqs. (2.17) and (2.25) is valid not only for Lagrangian, 
but also for general local functional C : 

R°T exp jC — T exp £[£] , (3.4) 

where 

£[£] = A UV °(T exp£-l). 
It is not difficult to see that £[£] is also a local functional: 

m = Y,f dqUC;q}] n (q). (3.5) 

n J 

In momentum representation its coefficient functions are formal series of the form: 
£n[£] q] = Yl ■ ■ ■ S / d qi-- - f dq N 5 (q + Y,q)j 

N=0 ni n N J J \ j=i / 

x P£ ) 1 ...n JV (9i • --9iv)An(5i) • • ->C niV (gjv), (3.6) 
where P are polynomials in g,. The UV poles are contained in the coefficients of the 
polynomials P. In the coordinate representation, x] are sums of local products of 
Ci(x) and their derivatives. 

The representation of £[£] in the form of (3.5) and (3.6) is in general not unique 
because a factor depending only on q as a whole can be included either into q] or 
into j n (q). But it can be made unique if one requires e.g. that V from (3.6) should not 
contain such factors. We assume that this condition is always fulfilled. 

If £ is a Lagrangian functional (cf. (3.3)) then such is £[£] as well, as can be seen 
from (3.6). 

So far we have done very little beyond rewriting the results of sect. 2 in a new form. 
But this exercise is by no means trivial because it exhibits the fact that the ^-mapping 
is universal in the sense that it depends only on the field content of the model but not 
on the interaction. We will see that the ^-mapping is a convenient tool for studying 
UV renormalization of interaction Lagrangians as well as composite operators. 



3.3 Inversion of the ^-mapping 

We are going to prove that there exists a mapping of the form of (3.5)-(3.6) 

with V replaced by some other polynomials V satisfying the same restrictions. More 
precisely, 

i\i-\L\\ = C and r'ffl^ (3-7) 

for any local functional C. The equations (3.7) allow one to write down the following 
inversion of (2.3): 

R°T exp r 1 [£] = T exp C. (3.8) 
This equation is one of the key results of the present paper. 



3 INVERSION OF THE R-OPERATION AND THE ^-MAPPING 



14 



3.4 The structure of the ^-mapping 



Let us study the structure of £[£]. (Throughout this section we denote A = A uv .) 
The terms that are linear in C can be extracted from £[£] as follows: 



i[C\ = A°(Te £ - 1) = A°(T£) + A°T (e c - C - l) 



(3.9) 



Then 



A°(T£) = AoT 



E / dq£ m (q)jm(q) = E/ dqC m (q)Ao(Tj m (q)) (3.10) 



and 



A°(Tj m (g)) = J2 Z m,njn(q)- 



(3.11) 



Indeed, j m (g) is the Fourier transform of a local products of fields so that the operation 
of T-product generates all possible tadpole graphs from j m (q). A replaces the loop 
integrals corresponding to such diagrams by counterterms, and the result is a linear 
combination of local operators, which is what (3.11) states. In any subtraction scheme 



= 1. 



(3.12) 



In the MS-scheme, z m>n for n ^ m are poles in D — 4 with numeric coefficients that are 
independent of the dimensional parameters. (Recall that we include masses into the 
local operators j.) 

Let us show that the matrix z m ^ n has an inverse. To this end we choose a special 
basis j n such that each j n is a monomial of fields and masses. Moreover, it is easy to 
see that the basis can be chosen to be ordered so that if the n-th operator is built of a 
lesser number of fields than the n'-th one (provided the dimensionalities in mass units 
and the numbers of full derivatives are the same) or if the dimensionality of the n-th 
operator is less than that of the n'-th one, then n < n' . In such a basis, m < n implies 
that z m . n = 0. Therefore, z mtn is a block-triangle matrix with units on the diagonal. 
It follows immediately that z~ l exists and has the same block-triangle structure as z, 
and the matrix elements of z~ x are polynomials of non-zero elements of z. 



3.5 Existence of £ 

Now we can solve the equation 

L = i\L\ (3.13) 

with respect to C. 

Rewrite (3.7) in the form of an equation for the sources: 



(3.14) 



3 INVERSION OF THE R-OPERATION AND THE ^-MAPPING 



15 



using the above information on the structure of £, extract the terms that are linear in 
£ on the r.h.s. of (3.14): 

£n(q) ='%2£m(q)Zm,n + 2€ n [£;q]- (3-15) 

m 

Multiplying (3.15) by z^ 1 , we can rewrite it as: 

£»(?) = 2 £m(q)z-] n - 51 26»[£; ?]^n" ( 3 - 16 ) 
m m 

Now the expression of C n in terms of C n can be obtained by iterating (3.16), taking as 
a starting value the first term on the r.h.s. of (3.16). k iterations allow one to obtain 
an expression of C n in terms of C n to order k in C n . Denoting the obtained solution of 
(3.16) as 

C = C 1 [£], (3-17) 
one can easily see that the mapping has the properties described in subsect. 3.3. 

3.6 A useful property of £ and 

If J is a source for the field cp then (2.32) can be rewritten as: 

Z[£ + <pJ]=Z[C]+<pJ. (3.18) 

It is not difficult to see that 

C 1 [C + <pJ]=C 1 [£]+<pJ- (3-19) 

3.7 ^-mapping and renormalization group 

It might be interesting to note that the above proof did not use any specific properties 
of the MS-scheme, though £ and do depend on the details of the i?-operation via 
V and C in (3.6). It follows that the above formalism offers an alternative way of 
studying renormalization group problems. 

For instance, let R a and H b be two i?-operations differing by the choice of subtrac- 
tion operators. Then: 

R a °T e c = T e uc] = R 6 °T e^K-Ml = R 6 oT e & [£] . (3.20) 

One sees that £ has the form of (3.5)-(3.6) with suitable polynomials V which — as 
can be easily understood — stay finite when the regularization is removed. Thus the 
structure of the renormalization group emerges very naturally 4 . 

4 It may be interesting to compare the above rather compact reasoning with the more cumbersome 
derivation of the renormalization group transformation given in [24] in the BPHZ framework. 
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The Zimmermann identities for the Green functions with local operator insertions 5 
can be easily obtained by performing suitable variations of (3.20) with respect to the 
classical sources C n entering C in (3.2) (cf. sect. 4 below). 

3.8 Practical calculations of R _1 . 

Let us derive a diagram- by-diagram recipe for calculating the mapping To this 
end it is convenient to introduce the inversion of the i?-operation, i.e. an operation 
R 1 such that its structure is the same as that of the i?-operation (see (2.1)) but 
the counterterms are evaluated according to a different rule in order to ensure that 
R-^R = 1 and R°R 1 = 1. 

By definition, for the inverse .R-operation one would have: 

R-^T e c = T exp £ _1 [£] . (3.21) 

There should exist the operator A~ x related to R _1 in the same way as A to R (see 
(2.1); recall that throughout this section A = A uv ). Then: 

r X [£] = A- 1 o(Te £ -l). (3.22) 

From the second equation in (3.7) the following condition on A -1 emerges: 

C = A _1 °[T exp A°(T e c - 1) - 1]. (3.23) 
Eq.(3.23) can be rewritten as: 

A- 1 o(Te £ - 1 - C) + (A-^C - C) 

= -A- x o[T expA°(Te £ - 1) -Te £ ], (3.24) 

where we have postulated linearity of A -1 . Note that the argument of A -1 on the 
r.h.s. of (3.24) is equal to (R — l)»(Te £ ). 

Assume that 

A-^[C] = C, (3.25) 

which means that A -1 is a unity operation on an elementary vertex. Then (3.24) can 
be rewritten as: 

A- 1 o(Te £ -l-£) = A- 1 o[(R-l)oTe £ ]. (3.26) 

Therefore, for (3.23) to hold, it is sufficient that the operator A~ x has the following 
properties: linearity; the property (3.25); equality to on unconnected and non-lPI 
diagrams; its action on 1PI diagrams should be defined by the following recursion: 

A _1 o[G] = -A- x o[R - l]oG. (3.27) 

5 It should be noted that we use the term "Zimmermann identities" in a slightly different sense than 
e.g. in [24]: we use it to denote the relations between the MS-renormalized and bare (unrenormalized) 
operators, without e-dependent factors like those considered in [24] for the purposes of studying 
anomalies. 
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The r.h.s. should be understood as follows. (R — 1) replaces G by a sum of integrals 
obtained from G by replacing one or more non-trivial UV-subgraphs by the correspond- 
ing UV-counterterm (see the definition of the i?-operation in subsect. 2.1). Therefore, 
A -1 on the r.h.s. acts on a linear combination of Feynman integrals with divergent co- 
efficients, all the integrals having one loop less as compared to G, while the divergent 
coefficients remain unaffected by A -1 because of its linearity. On one loop integrals: 

A -%] = -A-^Aofo] = -Ao[g], (3.28) 

and the recursion stops correctly. 

Eq.(3.27) provides a most convenient recipe for explicit calculations. For example, 
in the scalar gip A model one has: 

(3.29) 



and so on. 



4 Renormalization of multilocal operator insertions 

It is not difficult to generalize the above results to arbitrary Green functions of compos- 
ite operators. The main point here is that it is simply sufficient to consider T-exponents 
of arbitrary local functionals C as defined in (3.2) instead of Lagrangian functionals, 
and then derive expressions for specific composite (multi-) local operator insertions 
by performing suitable variations. In subsect. 4.1 we present some explicit formulae 
for the case when one wishes to express renormalized Green functions in terms of un- 
renormalized ones, and in subsect. 4.2 we consider the opposite case. Note that the 
case of one local operator insertion within the momentum subtraction scheme was first 
considered by Zimmermann [9]. 6 



6 see also the remark in the footnote in the last paragraph of subsect. 3.7. 
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4.1 Generalized Zimmermann identities 



Consider a Green function of arbitrary local operators (without loss of generality we 
take the local operators to belong to the basis introduced in subsect. 3.1): 



Ro < Tj 1 (g 1 )...j n (g n )e i > c 



(4.1; 



where jk are Fourier transformed local products of <p(x) and its derivatives as defined 
in (3.1). Using the formalism introduced in the preceding section, one gets: 



eq.(4.1) 



5 6 RTeA 

SCuApn) ' ' ' SCnM ) C=L: ^ =0 



C=L, ip=0 



5£ nN {p N ) ' ' ' <J£ m (pi) 



ML] 



(4.2) 



C=L 



V5=0 



The expression in the square brackets on the r.h.s. of (4.2) is a sum of products of 
terms like 



Tjn(p) • • .j in (p n )exp L 



(4.3) 



_5C in (p n ) ' ' ' <JAi(Pl) 

It is not difficult to see that (4.3) is a local functional which can be expanded in j as: 
eq.(4.3) = ^Jm(Pi + •••+*») 4"^;m(Pi ■ • -Pn), (4.4) 



where Z are polynomials in p and in the parameters of L. Z is independent of p if 
the dimensionality of j in units of mass is greater than that of the product . . . ji n . 
Moreover, since we include full derivatives into j, Z can not contain factors depending 
only on p 1 + . . . + p n \ in particular, Z^(p) are independent of p. 

Examples will help to understand the above. Denoting (4.3) as A(ji . . .j n ), one has 
(L R was defined in (2.26)): 



R° < T 3l e L > = (T A(ji) exp[L R ]) v=0 , 



(4.5) 



R° < T 3lj2 e L > = (T [A0'i)A0' 2 ) + A(j lj2 )] exp[L R ])^ =0 , 



(4.6) 



R° < Tjij2j 3 e L > = (T [A0'i)A(j 2 )A0' 3 ) + A(j lj2 )[A(j3) 



+ A(j 1 j3)A(j 2 ) + A(j 2 j 3 )A(j 1 ) + AUihJs)] exp[L R ]) v=0 . 



(4.7) 
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Since A(</?) = if and A(ipji ■ ■ ■ jk) =0 (this is because it is impossible to form a 
1PI graph with one "(^-insertion"), one can see that for j = <p the above formulae 
degenerate into 

Ro < T£( 9l ) . . .^{q n )e L > = (T£( 9l ) . . . y3(g n )e iR ^) v=0 , (4.8) 

which agrees with (2.33), as expected. 
Using (4.4), one can represent (4.5) as: 

Ro(T Up)z L ) = E Z^MUP) exp[L R ]) (4.9) 



(identities of this kind were first obtained by Zimmermann 

R°(T j 1 (p 1 )j 2 ( P 2)e L ) = E zJ 2 2 ) ;m (p 1 ,p 2 )(T ] m (pi + P2) exp[L R ]) 

m 

+ £ zg, E ^;i(T Upi)UP2) exp[L R ]), (4.10) 

etc. Note that all the sums in (4.9) and (4.10) are finite. Identities like (4.10) and 
their generalizations can be called (generalized) Zimmermann identities for multilocal 
operator insertions. 

4.2 Inverted Zimmermann identities 

The equations derived in the preceding subsection allow one to express renormal- 
ized Green functions in terms of unrenormalized ones. Using the inversion of the R- 
operation obtained in sect. 3, one can easily do the opposite, i.e. reexpress unrenorma- 
lized quantities in terms of renormalized ones. Thus, instead of (4.1)-(4.2) one has: 

< T~j 1 (q 1 )...] n (q n )e L > 

(5 5 

C=L, ip=0 



- - - Te £ 

5£ nN (p N ) <J£»i(pi) 



R°T 



8C nN (p N ) ' ' ' 5£ ni (pi) 

and instead of (4.9) one has: 



M , (4.11) 

C=L ) ip=0 



x(l) 

Tji(p)e L = £Z 1;m Ro(Tj m (p) exp T^D (4-12) 



~(i) 

(note that Z = [Z^]" 1 ), 



x(2) 

Tji(pi)i 2 (p 2 )e L = E^i,2;m(Pi,P2)R°(T j m (p! + p 2 ) exp C 1 ^]) 
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+ ^Z 1;m ^Z 1;n Ro(Tj m ( Pl )in(p 2 ) expT 1 ^]), (4-13) 

m n 
~(2) ~ (2) ~(1) 

(one can check that Z 12 . m (p 1 ,p 2 ) = - £j ^i,2;j(Pi,P2)^ ;m ). And so on. 

In a similar manner one can relate expressions renormalized in different renormal- 
ization schemes, cf. subsect. 3.7. 



ylS-OPERATION FOR MULTILOOP INTEGRALS. 
5 Heavy mass expansions 

We now turn to asymptotic expansions of Feynman diagrams. In this section we 
consider the case when the set of the heavy parameters contains only masses, which 
corresponds to studying low-energy effective Lagrangians to all orders in the inverse 
heavy mass. The reasoning in the more general case is very much the same, so that in 
the next section we will be able to concentrate only on specific complications due to 
large external momenta. 

In subsect. 5.1 the As-operation for unrenormalized (or UV-convergent) integrals is 
defined as an integrated version of the As-operation for products of singular functions 
introduced in [1], in subsect. 5.2 the combinatorial structure of the corresponding 
subgraphs is studied. In subsect. 5.3 the (^-functions corresponding to the counterterms 
for the IR-subgraphs are integrated out and in subsect. 5.4 the final form of the As- 
operation for 1PI integrals is presented. In subsect. 5.5 the obtained expression for 
the As-operation is extended to non-lPI integrals, in subsect. 5.6 its exponentiation 
is proved, and in subsect. 5.7 it is transformed to an explicitly convergent form using 
the inversion of the i?-operation. In subsects. 4.8 and 4.9 the results are extended 
to UV-renormalized diagrams and models, respectively. Explicit expressions for low- 
energy effective Lagrangians are presented to all orders of the inverse heavy mass. 
Subsect. 5.10 contains a simple example. 

5.1 ^-operation for integrals: definition 

Let us introduce notations and formulate the problem that we are going to consider 
here. Our notations on the whole will be consistent with subsect. 3.1 of [1]. 

Let T be an /-loop 1PI Feynman diagram, and let p = (pi . ■ .pi) denote its inte- 
gration momenta. We assume in this section that V depends on the light and heavy 
masses m and M, and that all the external momenta of T are light (and denoted as k). 
All M are non-zero and independent of the expansion parameter denoted as k, while 
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m,k = 0(n) and some masses m may be equal to zero. Unlike the abstract analysis of 
[1], here it will be convenient to show the dependence on m, k and M explicitly. 

Thus, the unrenormalized integrand is denoted as T(p, m, k, M). We assume that Y 
is UV- convergent, i.e. T(p, . . .) is absolutely integrable over p in infinite limits. Denote 
the integrated diagram as: 

T(k,m,M)= J dpT(p,k,m,M). (5.1) 

In [1] we introduced the As-operation for products of singular functions like T(p, k, m, M) 
which yields asymptotic expansions of such products with respect to the light parame- 
ters in the sense of the distribution theory. It is natural to expect (for more details see 
subsect. 5.8; a fuller justification is given in [17]) that if (5.1) is UV-convergent, then 
the asymptotic expansion of T(k,m, M) can be obtained by expanding its integrand 
r(p, k,m, M) using the As-operation of [1] and then performing termwise integration 
of the resulting series over p. 

Let us introduce an operation that acts on integrated Feynman diagrams in accor- 
dance with the above recipe. This new operation will be referred to as the As -operation 
for integrals. To distinguish the two types of As-operations, one could use e.g. the 
notations As pro( j and As int . But as it is normally clear whether the As-operation acts 
on a product of singular functions or on an integral, we will use the same notation As 
without subscripts in both cases. Formally, one has: 

As°T(k,m,M)= J dpAs°T(p, k, m, M). (5.2) 

Let us exhibit the structure of the As-operation on the r.h.s. of (5.2). The general 
expression for it was given in [1]. Specifying it to the case considered here, one obtains: 

As°r(j9, k,m, M) 

= T m , fc °rGo, k, m, M) + £[Eo 7 (p, k, m)] x [T m . fe or\ 7 (p, k, m, M)], (5.3) 

7 

where the notations used are as follows. 

The operation T m ,k performs the Taylor expansion in powers of the light parame- 
ters. Summation runs over all IR-subgraphs of Y. Specifying to the present case the 
definitions of subsects. 9.2-3 of [1], one sees that an IR-subgraph 7 can be described 
as follows: (i) the set of lines of 7 is a full subset of the light lines of T, i.e. if one puts 
to zero the momenta flowing through all the lines of 7 together with the light external 
momenta k, then no light line which does not belong to 7 should have the momentum 
flowing through it nullified owing to momentum conservation (cf. [1]); (ii) any vertex 
of T whose incident lines all belong to 7, is included into 7. The product of the lines 
and vertices of 7 is independent of M and is denoted as j(p, k, m). 

The product T\7 in (5.3) comprises all the factors from V that do not belong to 7. 
Graphically, T\7 is obtained from V by deleting all the lines and vertices of 7. 
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The operation E (expressions for it were derived in [1]) applied to an IR subgraph 
returns a linear combination of counterterms that are proportional to (^-functions of p. 
Our first aim here is to integrate out those 5-functions explicitly. 

5.2 Structure of T\ 7 . 

It is not difficult to verify that in the case under consideration (no heavy external 
momenta) the IR-subgraphs are uniquely characterized by the properties of their com- 
plements. Indeed, denote the connected components of T\7 as hi. These subgraphs 
(which will be referred to as heavy knots) have the following properties: 

(i) since 7 consists of only light lines, each hi has only light external lines while all 
the heavy lines of T are hidden within all hf, 

(ii) since 7 is full, none of hi can have a light line whose momentum vanishes due 
to momentum conservation when one puts to zero all the external momenta of hi. This 
is an analytical formulation. Graphically, it implies that each hi is "1PI with respect 
to light lines" which means that it cannot be divided into two disconnected graphs 
(which may consist of a single vertex) by cutting any one of its light lines. 

Note that a heavy knot hi may consist of a single heavy line (including both vertices 
to which it is attached). 

One has: 

r\7 = n^- ( 5 - 4 ) 

i 

In the final form of the As-operation for integrals, the heavy knots will play the same 
role as UV-subgraphs in the expression (2.1) for the i?-operation. 

Example. Consider the diagram V as shown in Fig. la. Fat lines correspond to 
propagators of heavy particles, the rest to light particles. In Figs. lb-Id various IR 
subgraphs are shown with dashed lines. Heavy knots in each case are clearly visible as 
connected components built of solid — normal and fat — lines. 

5.3 Integrating ^-functions 

Before substituting explicit expressions for E°7 into (5.3), the dependence of each 
term in (5.3) on p should be studied. One can see that 7 does not depend on those 
components of p that correspond to the loops of T\7. Let us denote the collection of 
such momenta as Pryy,mt- Furthermore, each loop of T\7 belongs to only one heavy 
knot hi from the decomposition (5.4), so that: 

Pr\ 7 ,int = {Pi,mt}i- ( 5 - 5 ) 

The remaining momentum components of p form exactly the proper variables of 7 
introduced in [1] and denoted as p 7 , so that 

P= (p 7 >Pr\7,int)- (5-6) 
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Note that T\7 is independent of those components of p 7 that correspond to the loops 
of 7. Denote as Pryy,ext those components of 7 on which T\7 does depend. 

Substituting (5.3) into (5.2) and using the explicit expressions for E°7 derived in 
[1], we can rewrite (5.2) as: 

As°r(£;, m,M) = T m , k °T(k, m, M) 



7 a 



x T 



m,k 



r \7(Pr\7,cxt,Pr\7,int, k, to, M)J }, (5.7) 

The polynomials V a and the ^-functions 5 a in (5.7) form full dual sets (cf. [1]) so that 
the sum over a represents the Taylor expansion: 



J2 v ~fAPy) J dp y 5^ a (p 7 ) x /(p 7 ) = T pL( of(p^). 



(5.8) 



(It does not matter that the expression that plays the role of / in (5.7) depends not 
on all p 7 but only on some of the components of p 1 , namely, on pryy.ext-) Using (5.8) 
and renaming p^ cxt -> p 7)6Xt , we get: 

As°r(/c, m, M) = T mr T(k, m, M) 

(5.9) 



+ H { j d Pj(Pv k ' m ) X T ™,fe° T Pr\ 7 ,cxt° [ r \7(pr\7,ext,Pr\7,int, k, 771, M) 

Now recall the representation of T\7 in terms of heavy knots (5.4), (5.5). It is easy to 
see that the last line in (5.9) can be rewritten as 



II T ° h i (PiM ' Pi,«t ,k,m,M), 



(5.10) 



where Pi, ex t comprises those components of pr\7,cxt on which the heavy knot hi depends. 
The operation r in (5.10) Taylor-expands in m and in all the momenta that are external 
with respect to the subgraph hi, i.e. both p^ext and the corresponding components of k. 
It should be stressed that each factor in (5.10) is independent of the rest of the original 
diagram Y. Therefore, we arrive at the following expression of the As-operation for 
integrated Feynman diagrams: 



(5.11) 



A S or = Tor+]T x [r/IW- 

{hi} V i 1 i 

The meaning of this formula is as follows. One enumerates all sets of pairwise non- 
intersecting heavy knots hi and contracts each heavy knot hi to a point, replacing it 
by a formal series in its external momenta which is obtained by Taylor-expanding the 
corresponding Feynman integral. 

To avoid confusion, it should be understood that the IR-counterterms E°7 in (5.3) 
have transformed into the expression in the square brackets in (5.11), while the "coun- 
terterms" r°hi in (5.11) have originated from the expression T\7 in (5.3). In other 
words, "counterterms" and "non-counterterms" have changed places in the transition 
from (5.3) to (5.11). 
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5.4 ^-operation for integrals: the final form 

Eq.(5.11) can be rewritten in the following form which is analogous to the representa- 
tion (2.1) of the i?-operation: 



The summation here runs over all sets of pairwise non-intersecting subgraphs hi of T 
such that each vertex of T belongs to one of hi. (A subgraph h is defined as a subset 
of vertices of T together with some lines of T which are connected by both ends to 
the vertices included into h. Subgraphs are non-intersecting if they have no common 
vertices and, consequently, no common lines.) One of such sets consists of a single 
subgraph h — T which corresponds to the first term r°T on the r.h.s. of (5.11). The 
operation A as coincides with r on heavy knots, is a unity operation on single vertices 
without loops, and returns zero otherwise. Setting M = oo is a formal expression of 
the fact that one should discard from the r.h.s. those terms in which some heavy lines 
remained outside all hi so that (5.11) could be restored from (5.12). 

There is no need to introduce special rules for heavy knots containing tadpoles 
made of light lines (cf. Fig. 2), our formalism takes them into account correctly. For 
practical purposes it is sufficient to note that such heavy knots are nullified by the 
operation r due to the well-known property of the dimensional regular izat ion to nullify 
momentum integrals without dimensional parameters. 

The combinatorial resemblance of (5.12) and the expression (2.1) for the i?-operation 
is nearly complete except for the minor difference described above. However, this dif- 
ference can be easily taken into account to make the reasoning of sects. 1-3 to be 
applicable to the As-operation. But first expansions of diagrams T other than 1PI 
should be considered. 

5.5 ^s-operation for non-lPI diagrams 

We have derived explicit expressions for the operation As on 1PI integrals. Remarkably, 
one can expand arbitrary — not necessarily 1PI — integrals using the same expression 



For the operation As to be applicable to any graph T, its structure should satisfy 
certain a priori requirements. Let us first enumerate them. 

If T is disconnected, i.e. there are two subgraphs Ti and r 2 such that r = T\ x r 2 , 
then the expansion of T is a product of the expansions of Ti and T 2 , i.e. 




(5.12) 



(5.12). 



As°(r! x T 2 ) = (As°Ti) x (As°r 2 ). 



(5.13) 



If r falls into two parts connected by only one light line (the corresponding propa- 
gator is denoted as A), i-e. T — T 1 x Di x T 2 , then 

Aso(r! x a x r 2 ) = (As<.ri) x a x (Asor 2 ). (5.14) 
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This is due to the fact that both the mass of this line and the momentum flowing 
through it (which is a combination of the light external momenta of T) are expan- 
sion parameters so that Di is proportional to a negative power of k and need not be 
expanded. 

If T consists of two subgraphs connected by one heavy line (the corresponding 
propagator is denoted as D h ), i.e. r = r\ x D h x T 2 , then 

Aso(i\ x D h x r 2 ) = (AscrO x (T k oD h ) x (As°r 2 ). (5.15) 

(Recall that T fc Taylor-expands in k, and that the momentum flowing through D h is a 
combination of the light external momenta of V.) 

Now, note that the definition (5.12) does not use lPI-ness of V . Therefore, let us 
define the operation As by (5.12) on any graph Y. Let us prove that the operation thus 
defined satisfies the above a priori criteria (5.13)-(5.15). 

Eq.(5.13) is true because each hi from (5.12) lies strictly within either Ti or T 2 , so 
that each partition of Y consists of a partition of Y\ and a partition of T 2 , and a sum 
over all partitions of Y is equivalent to a sum over partitions of Ti times a sum over 
partitions of T 2 . 

The same argument proves (5.14) because Di cannot belong to any of hi without 
making it 1PI with respect to light lines, which implies that each hi with non-zero A as 
lies strictly within either Ti or T 2 . 

To prove (5.15) we note that Dh must always belong to one of the hi. Denote this 
heavy knot as h. Then h = h' x Dh x h" (where h' and h" each consist of at least one 
vertex) and 

A as o(ti xD h x h") = (A as oti) x (T k oD h ) x (A as oh"). (5.16) 

Then one applies the same argument as above. 

An induction with respect to the number of 1PI components of Y enables one to 
conclude that (5.12) remains correct on all Feynman diagrams. 

5.6 Exponentiation of A as . 

Now we are in a position to apply the combinatorial results of sects. 1-3 to (5.12). 

Let ip and $ be the fields of the light and heavy particles that have the masses 
m and M, respectively. Let L = L((p, <3>) be an interaction Lagrangian. Consider the 
generating functional of Green functions of the light fields < Texp[L(</?, $) + tpJ] >o • 
It is assumed that all the external momenta of the Green functions are light, i.e. O(k). 
We ignore the UV divergences for a while, because we are now interested only in the 
formal global structure of the operation As; one can simply assume that there are no 
UV divergences at all in the model described by L. 

As was noted above, the combinatorial structures of the As- and i?-operations are 
very similar. The similarity goes far enough to allow one to write down the following 
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analogue of (2.31): 



As° < Texp[L(<^, $) + ipj] > = < Texp[L eff)barc (^) + tpj] >„, 



(5.17) 



where 



L eff ,barc(^) = A as ° W e L{ ^ - 1 



J$=0 



(5.18) 



Setting $ to zero in (5.18) is formally equivalent to setting M = oo in (5.12) and is an 
expression of the fact that each heavy line in (5.11) must belong to one of the heavy 
knots. 

The algorithm of (5.18) is as follows: 

(i) all the diagrams contributing to T exp L(tp, <3>) are evaluated; 

(ii) A as nullifies all the terms except those which correspond to connected diagrams 
that are 1PI with respect to light lines; 

(in) on single vertices, A^ acts as a unit operation, i.e. 



(iv) the remaining terms are exactly the heavy knots and can be represented as 
(2.28), where FI also depends on m and M. A as performs the formal Taylor expansion 
of FI with respect to m and all the external momenta of such terms, i.e. pi in (2.28). 

Note that the recipe for A as allows a remarkably "fool-proof reformulation (which 
may be useful e.g. in writing computer programs): it is sufficient to apply the r- 
operation to all connected terms generated by the T-exponent. Then the terms that 
should be discarded will either develop singularities 1/0, or be nullified owing to the 
property of the dimensional regularization to nullify integrals without dimensional pa- 
rameters like 



It is no less remarkable that this "fool-proof reformulation stays valid in the most 
general case comprising heavy external momenta with arbitrary linear restrictions. 

5.7 Explicitly convergent form for the heavy mass expansion 

As was noted in [1], our method of expansion leads to expressions containing spurious 
divergences in different terms which cancel in the final result. Let us reorganize our 
asymptotic expansions in such a way as to reexpress them in an explicitly convergent 
form. In (5.18), L e ff,barc(y?) is a Lagrangian functional (in terms of subsect. 3.1) of 
the field (p whose coefficients ( "coupling constants" ) are functions of the heavy masses. 
Irrespective of the specific values of these functions, the diagrams on the r.h.s. of (5.17) 
may contain UV divergences. Using (3.11), one can explicitly extract the -R-operation 
on the r.h.s. of (5.17) and rewrite it as: 



A as o[L(v?,<l»)]$ = o 



[L((p, $)]d> =0 ; 



(5.19) 




As° < Texp[L(<£, $) + cpj] > = R° < Texp[L cff (y?) + tpj] > , (5.20) 
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where 



Le&iv) = £ ^eff.bare (<£>)]■ 



(5.21) 



Since the mapping £ 1 transforms a Lagrangian functional again into a Lagrangian 
functional, L e ff(y?) can be represented as 



where j n (x) are defined in subsect. 3.1 and g denotes the set of coupling constants of 
the original interaction Lagrangian L; the "effective couplings" g e s,n have the form: 



where the summation over powers of g reminds one that we are working within pertur- 
bative framework; g e s,n,i(t) have a pole singularity at e ~ of a finite order depending 
on n and I, while the integer exponents d n (some which for a few lowest values of n 
may be negative) are determined on dimensional grounds. 

Indeed, the coefficients g e ff,n are polynomial expressions built of UV counterterms 
(each of which is just a polynomial in 1/e) and heavy knots which are Taylor-expanded 
in powers of their external momenta and m. The resulting loop integrals depend only 
on M, and on dimensional grounds this dependence has the form M d ~ 2d x c(e) where 
d! is an integer and /' the corresponding loop number, while c(e) has a pole singularity 
l/e l (the latter information is strictly speaking not necessary). The powers of momenta 
and the light masses m resulting from the above Taylor expansion are included into j 
(recall the definition in subsect. 3.1). Recalling that the UV renormalization parameter 
\i enters all expressions via integer powers of /i 2e , one arrives at (5.22) and (5.23). Note 
that the range of summation over / is finite for each n and k. 

Now one can easily convince oneself that the parameters of the effective Lagrangian 
are finite at e = 0. Indeed, the initial expression (the l.h.s. of (5.20)) is finite at e = by 
construction in each order of the perturbation theory (in the absence of UV divergences 
in the initial model — recall the assumption at the beginning of this section). When 
transformed into the form of the r.h.s. of (5.20), it is, of course, still finite but has now 
the form of UV-renormalized (and therefore, finite) loop integrals times products of 
9eS,n{g, M, /i, e), in each order of the perturbation theory. Taking into account linear 
independence of j n , one concludes that all g e s, n (g, M, /i, e) are finite at e = 0. 

5.8 ^-operation on UV-divergent diagrams 

Eq.(5.2) was used to define the As-operation on UV-convergent diagrams. However, 
that definition can be formally extended to the UV-divergent case (then, of course, 
the UV divergences should be regularized; but the dimensional regularization which 
is used at the intermediate steps of the As-operation for products as derived in [1], 




n 



(5.22) 




(5.23) 
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regularizes UV divergences as well; therefore, no new regulators are really needed). 
The combinatorial results (5.17)-(5.21) remain valid though now the effective couplings 
9eS,n(g, M, fi, e) are not finite at e ~ because they inherit the UV divergences of the 
unrenormalized Green functions on the l.h.s. of (5.17). 

Let R°r(/c, m, M) be a renormalized graph (the dependence on the renormalization 
parameter \i is implicit; note that we are dealing here with integrated graphs, so that 
there is no dependence on the loop momenta). The i?-operation was defined by (2.1), 
which can be represented as follows: 



where % enumerates all the individual terms in the sum in (2.1), Z^ are products of 
divergent coefficients of the corresponding UV counterterms and are polynomials of all 
their arguments, Y { are the (unrenormalized) diagrams obtained from T by shrinking 
the corresponding UV-subgraphs to points. 

Naively, to expand (5.24) one would expand each term on the r.h.s.: the factors 
Z i are polynomials in m and are already "expanded", and Tj can be expanded by a 
straightforward application of the definition (5.2). Indeed, all the intricacies of the 
As-operation are essentially aimed at extracting the terms that are non-analytical in 
light parameters, and the very fact of the polynomial dependence of UV counterterms 
on masses and momenta makes one expect that UV renormalization in the MS-scheme 
and the yls-expansion are, in a sense, "orthogonal" (cf. the analysis of [17]). 

Motivated by the above, we define the As-operation on renormalized graphs by 
applying the As-operation (5.2) termwise to Ti on the r.h.s. of (5.24) as follows: 



We assert that the expansion thus obtained is a correct asymptotic expansion for 
R°r(/c, m,M). This is not obvious, but for that matter neither is obvious the UV 
finiteness of the r.h.s. of (5.24) or even (2.1). 

For a practically oriented reader it may be sufficient to learn the following. On the 
one hand, the final expressions for the coefficient functions of short-distance expansions 
(see subsect. 6.3) imply highly non-trivial cancellations of UV and IR divergences 
between different terms which cannot be understood without the underlying theory. 
Such cancellations were checked by explicit three- loop calculations in [11] and in several 
other calculations (see e.g. [18]). On the other hand, from the point of view of our 
methods there is no essential difference between the short-distance operator-product 
expansion and the most general Euclidean regimes with heavy masses, and all the non- 
trivial patterns of interaction of UV and IR divergences manifest themselves at the 
three-loop level. One should also take into account that we are dealing with relatively 
primitive although increasingly cumbersome integrals of rational functions where all 
the effects are governed by the power counting (this point is emphasized and exhibited 
in the regularization independent formalism developed in the parallel series of papers 



K°r(k,m,M) = ^Zi(l/e,m,M)Yi(k,m, M), 



(5.24) 



AsoRor(A;, m, M) = ]T Z^l/e, m, M)As°Ti(k, m, M). 



(5.25) 
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[14]-IV), so that there is no room for pathologies. Therefore, the explicit calculations 
of [11] and [18] provide a sufficient evidence in favour of the above assumption. 

For a formally oriented reader we note that as was proved in [15], the i?-operation 
can be represented in such a form that its UV finiteness becomes obvious. Moreover, 
the reasoning of [17] demonstrates that the proof of correctness of the above defini- 
tion reduces to verifying commutativity of two As-operations — one corresponding to 
the expansion under consideration, the other being implicit in the construction of UV 
subtractions along the lines of [15]. Such commutativity is essentially due to commuta- 
tivity of the corresponding formal expansions of the individual factors, and verification 
of the fact that it is indeed inherited by As-operations on the entire integrands is a 
rather straightforward technical excercise. An informal discussion of the proofs of [17] 



In what follows, we will concentrate on the combinatorial structure of the definition 



(Note that the aggregate operation As°R defined in (5.25) was denoted as As and 
studied as a whole in our original publication [2], and the starting point there was the 
so-called EA-expansion. Such a way of reasoning resulted in superfluous complications 
which are avoided in the present exposition.) 

5.9 v4s-operation in models with UV divergences 

If the model described by the interaction Lagrangian L(ip, $) possesses UV divergences, 
then the UV-renormalized generating functional of Green functions can be represented 



with Lr(<p, $) defined analogously to (2.33). 

Since the As-operation is defined to act formally on each term of the renormalized 
diagrams, one can apply it to both sides of (5.26). Then one can simply use the results 
obtained for the case of the models without UV divergences (see (5. 17)— (5.21)) . Our 
final formulae are as follows: 



R° < Texp[L(<^, $) + tpj] > ~ As°R° < Texp[L(y?, $) + tpj] > 

M^oo 

= RoR-^AsoRo < Texp[L(v?, $) + ipJ] > = R° < T exp[L cff (v?) + ipJ] > , (5.27) 



is [16]. 



(5.25). 



as (cf. (2.33)): 



Ro < Texp[L(y?, $) + yj] > =< Texp[L fl (<p, $) + <pj] > . 



(5.26) 



where 



A„v» Te Lc «- barc(v) - 1 , 



(5.28) 



-^cff,barc(V 9 ) 



L R {p, $) 




(5.29) 



(5.30) 
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The algorithms of the three A-operations have been described in detail in sects. 1.1, 
2.8, and 4.4. 

Eq.(5.27) means that the effects of virtual presence of heavy particles on the effective 
low energy theory of light particles can be described by an effective Lagrangian L e g 
(explicitly given by (5.28)-(5.30)) to all orders in 1/M. Our derivation of this result 
is valid within the MS-like renormalization schemes and in all models including those 
with massless particles like QCD. 

Note that the i?-operations in the original model and in the effective model need 
not be the same; this can be equivalently described by saying that the renormalization 
parameters in the two models may be different. In such a case the parameters of L e g 
will depend on both renormalization parameters. 

It should also be recalled that the As-operation derived in [1] is essentially unique 
because it leads to series in powers and logarithms of the expansion parameters. This 
uniqueness property (up to the choice of renormalization schemes) is inherited by 
(5.27)-(5.30), whatever the guise in which such expansions might appear when de- 
rived by alternative methods. Our form of presentation is dictated by computational 
convenience. 

The property of uniqueness also facilitates the study of gauge properties of L c s (cf. 
[10]). 



5.10 Example 



To get a feeling of how the above formulae work, consider QED with a light electron ip 
and a heavy muon The well-known interaction Lagrangian is denoted as L(ip, A) 
where A is the photon field. Restricting ourselves to the two-photon sector in the 
1-loop approximation, we get (fat lines correspond to the muon): 

L R (ip, A) = A uv ° [T e^DOW) _ x 



and, finally, 



= A uv ° 



/ 



5Z i / dx 



+ A uv c 
--FL(x) 



4~ /xfV"vj + • • • 5 

^bare,eff(^,A)=A as o[Te LR ^ A )-l 



L R (ip, A) + A as o 



+ 



</>=o 
+ ..., 



L eS (ip, A) = [T e Lcff - baroW, ' A) - 1 

= L eff)b are(^, -4) + Ayy 



(5.31) 



(5.32) 



(5.33) 
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Since 



we obtain: 



L eS {ij,A) = A a 



-A uv ° 



+ A uv ° 



(5.34) 



+ 



= 5z{M 2 ,fi 2 )i J dx 
where the dots denote other field structures, and 



- 4 F » 



+ ..., 



Sz(M 2 ,/j, 2 ' 



16tt 2 



Cl + c 2 log(M7//) + 0(1/M 2 ) 



(5.35) 



(5.36) 



Note that cancellation of the UV divergences in the final result is obvious even without 
explicit calculations. 



6 Generalized operator-product expansions 

In this section we turn to the case when the set of heavy parameters includes external 
momenta. Recall that we regard the diagrams and Green functions to be expanded as 
distributions with respect to the heavy external momenta [1]. 

We start by considering in subsects. 6.1-6.3 the case corresponding to the usual 
short-distance operator-product expansion, i.e. when there are no heavy masses and 
no linear restrictions on the heavy external momenta except the overall momentum 
conservation. In subsect. 6.1 notations are introduced and the As-operation for inte- 
grals is derived for the case under study, in subsect. 6.2 the global structure of the 
As-operation as applied to Green functions is studied, and in subsect. 6.3 the final 
results are presented. Modifications due to heavy masses are studied in subsect. 6.4. 

In subsect. 6.5 we impose the so-called "natural" linear restrictions on the heavy 
momenta, and in subsect. 6.6 the corresponding version of the As-operation and the 
expansions for Green functions are presented. In subsect. 6.7 the structure of the 
general formulae is explained with a simple example. 

In subsect. 6.8 we briefly discuss the contact terms in the obtained expansions. 
The section is concluded in subsect. 6.9 with a discussion of the most general linear 
restrictions imposed on heavy momenta; it is pointed out that the contact terms in 
such a case may contain the so-called paralocal operators. 



6.1 ^-operation and heavy momenta 



Let us assume that all the particles of the model are light, but some of the external 
momenta of the Green function to be expanded are heavy. 
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More precisely, let j numerate the heavy external lines of the diagram T, then Qj 
are the corresponding ingoing external momenta. There are also some light external 
momenta collectively denoted as k, which are of order O(k) together with the masses 
m (recall that k is our standard notation of the expansion parameter [1]). The fact 
that Q are heavy is formally expressed as 

Q = 0(1) as k^O. (6.1) 

One cannot assume all Qj to be independent of k because of the momentum conserva- 
tion which reads: 

£Q = -E* = °(*)- (6.2) 

Therefore, let us introduce the "heavy" components of Q that are independent of k : 

Qj = Qj\ K =o- (6.3) 

Q can be represented as: 

Qj = Qj + qj, (6.4) 

where qj = O(k) are linear combinations of k. 

The momentum conservation should hold separately for heavy and light compo- 
nents, so that: 

E& = 0. (6-5) 

3 

We have assumed that there are no other restrictions on Q except (6.5). Considering 
the diagram T as a distribution with respect to the heavy momenta means that we 
expand expressions of the form 

T{k,m,F)= J dQ' F(Q) J dpT(p, Q + q, k,m), (6.6) 

where F is a smooth test function independent of k, and the integration runs over the 
manifold described by (6.5). 

It is convenient to represent F as a vertex attached to the vertices corresponding 
to Q by the lines that are heavy by definition — see Fig. 3. 

To expand T(k, m, F), one repeats the reasoning of subsects. 4.1-4.4 and arrives at 
the following equation instead of (5.12): 

Asor = £(A as °/>) x [Y/h]. (6.7) 

h 

Now there is only one heavy knot shrunk to the point in each term in the sum. Its 
description — owing to the above agreement that the lines connecting the F-vertex with 
T are heavy — coincides with that given in subsect. 5.2: it must be 1PI with respect to 
light lines. But it may be easier to follow the "fool-proof recipe given at the end of 
subsect. 5.6, which remains valid here. 

A reasoning similar to that of subsect. 5.5 shows that the above formula is also 
valid in the case of a non-lPI graph T; in particular, it is valid for disconnected graphs. 
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6.2 ^-operation on Green functions 

Let us now turn to Green functions. Let Hj{x) be local products of light fields (p(x) 
and their derivatives. Consider the following generating functional of Green functions: 

<T{n^(g j )e £ }>o, (6.8) 

3 

where the tilde marks Fourier transforms and £ is a local functional (see the definition 
in subsect. 3.1). To obtain specific correlators from (6.8), it is sufficient to perform 
suitable variations with respect to the coefficient functions of C and replace C by the 
Lagrangian L. The momenta corresponding to any additional operator insertions are 
light, i.e. O(k), by definition. 

Without loss of generality, the test function is introduced as follows: 

G(F,C) =< T{F*He c } >„, (6.9) 

where 

F*H = J dQ f F(Q)Y[H j (Q j + qj), (6.10) 

3 

where the momenta Qj and qj are the same as defined in subsect. 6.1. 

The sum of all qj should in general not be taken to be zero in (6.9) in order to 
get rid of the disconnected diagrams contributing to (6.9) which normally are of no 
interest in phenomenological applications. If, however, the kinematics of the problem 
require that J2<lj — for the connected component, one can take the corresponding 
limit termwise in the final expansion. Taking such limits commutes with the expansion 
procedure unless there are connected diagrams like the one shown in Fig. 4 where the 
wavy lines correspond to massless particles and may give rise to an infrared diver- 
gence. Whether or not the problem under study allows such terms can be checked by a 
straightforward inspection. However, it is still possible to regularize such contributions 
e.g. by introducing a non-zero mass (which, of course, should be considered as light in 
the expansion procedure); the Q-dependent part of the final expansion as obtained by 
our methods will be insensitive to the light mass structure of the model (see below). 

Now, using (6.7) and reasoning as in subsect. 5.6, we obtain the following analogue 
of (5.17): 

As°G(F, C) =< T {A as o [TF * He c ] e c ) > . (6.11) 

The algorithm for evaluating A as [TF * HexpC] is exactly the same as described in 
subsect. 5.6 (recall that according to the graphical conventions introduced after (6.6), 
the lines connecting the F-vertex with the proper Feynman diagram are heavy by 
definition; note also that the 5-function which expresses momentum conservation in 
the relevant diagrams contains a sum of only light momenta due to (6.5), so that 
a formal expansion in light parameters will not affect it; if one wishes, one could 
introduce a formal integration in some of the qj to get rid of it completely). One sees 
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that A as [TF * HexpC] is a local functional in the sense of subsect. 3.1, so that (cf. 
subsect. 6.8): 

A as o [TF * He c ] = ]T C h Z c AF, C)j n (q), (6.12) 

n 

where q = J2Qj an d 

C*£(F,C)= J dQ'F(Q)C^(Q,C). (6.13) 

We have introduced the subscript "bare" to indicate that the C's contain divergences, 
so that a procedure analogous to what was described in subsect. 5.7 is needed, in order 
to transform (6.11) to an explicitly convergent form. 

To accomplish this, one proceeds as follows: 

As°G(F, C)=Y <%g(F, C)<T {] n (q)e C } >o 

n 

= E^(F,£)^R.<T^-M> 
= E ill Ch"(F, £)Zn,Jl R- < T {j,„(g)eS-'M} >„, (6.14) 

m \ n ) 

where we have used the fact that defined in sect. 3 is a local functional. 

To complete the reasoning, one should take into account that the initial expression 
on which the As-operation acts may contain UV divergences and be renormalized via 
the i?-operation, as follows: 

G R (F,C) = Ro < T{F*He c } > . (6.15) 

The arguments justifying correctness of a straightforward application of As to renor- 
malized expressions are the same as in subsect. 5.8. 



6.3 Short-distance expansions 

Now let us exhibit the general structure of the expressions for the case when C = L+ipJ 
where L is a Lagrangian functional. 

Using the results of sect. 4 (cf. (4.9)), one represents (6.15) as: 

Gr(F, L + V J) = Y<^{f* Z n , a * JJ Lr+ ^} > , (6.16) 

a 

where J a are multilocal products of the currents j n defined in subsect. 3.1, Z is a "ma- 
trix" of divergent coefficients, and the *'s denote contractions over the heavy momenta. 
Lr is defined in (2.26). Applying As to both sides and using (6.14), one gets: 

As°G R (F,L + V J)= I dQ'F(Q)YcUQ,g,v)R°<T{] m ((l)e [L+vJ] } >o, (6.17) 
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where 

/ dQ'F(Q)c m (Q,9^) = J2J2 C jZn( F * z n,a,L R )Z^ m (6.18) 

J an 

are linear functionals of the test function F that are finite in the limit e — > (which 
is proved similarly to the reasoning at the end of subsect. 5.7); g are the coupling 
constants of the Lagrangian L, and ji is the renormalization parameter. The expansion 
(6.17) has the form of a familiar operator product expansion at short distances. 

It should be stressed that (6.18) — despite its cumbersome appearance — describes 
a fully constructive algorithm for getting explicit expressions for the coefficient func- 
tions C m (Q, g, fi). This algorithm is equivalent to the one described in [11], where an 
assumption on existence and properties of the operator-product expansion in the MS- 
scheme was made. Namely, in [11] it was assumed that the coefficient functions of such 
OPE are analytical in masses. One can easily see that the expressions (6.17)-(6.18) 
prove that assumption: indeed, the coefficient functions as defined above have turned 
out to be independent of to, but the local operators j have been allowed to contain 
non-negative integer powers of to; if the local operators are built of only fields, then 
the coefficient functions would become analytical in masses. Non-trivial examples of 
two- and three-loop calculations of coefficient functions using this algorithm together 
with explicit formulae can be found in [11] and [18]. 

Finally, note that if one wishes to consider Green functions of composite operators, 
then suitable variations in C should be introduced into (6.14), (6.15) etc., giving rise to 
additional terms in the final result. Thus, if one considers the expansion of a correlator 
of the form 

< T {j n (Q)j m (-Q + q)]i(k)e L } > , (6.19) 
then the final expansion will contain the correlators 

<T{j a (g)j b (A;) e L }>o and < T {j a (q + k)e L ) > . (6.20) 

The latter expression is a typical "vacuum condensate" , i.e. a vacuum expectation value 
of a local operator (cf. [12]). It is not equal to zero, even within perturbation theory, 
because the normal ordering of local products of fields is not used in the MS-scheme. 



6.4 Effects of heavy particles on operator expansion 



If there are heavy masses in the model, then the above results will get modified in the 
following way. The expression (5.12) for the As-operation will remain valid provided 
one includes the heavy knots corresponding to the F-vertex that were described after 
(6.7). Note that the "fool-proof recipe of subsect. 5.2 is still valid here. Instead of 
(6.11) one will have: 



As°G(F,£) =< T{A a 



TF * He L 



exp A as ° 



Te L -1 



}>« 



(6.21) 



And the final result (6.17) will take the form 



As°G R (F,L(^, $) + y?J) 
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= / dQ'F(Q)J2c m (Q,M,g,^q)xRo < Ti i j m (J2q)exp[L cfi (ip) + V J}}>o, 

(6.22) 

where L ef[ ({p) is defined in (5.21), while c m are analytical in q. Note that the operators 
j m are here built of the light fields only, while the currents entering into the multilocal 
operator H in the initial expression are allowed to contain heavy fields as well. 



6.5 Natural restrictions on heavy momenta 

Let us turn to the case when there are linear restrictions imposed on the heavy external 
momenta other than the overall momentum conservation. Of immediate phenomeno- 
logical importance are the so-called natural restrictions [3], [1]. In terms of position 
space, one can describe them as follows: all the "heavy" operator insertions are ar- 
ranged into several groups, the distances within each group tend to zero while the 
distances between groups stay finite (cf. Fig. 5). More precise definitions in terms of 
momentum representation are presented below. 

Consider (6.8) and let the "heavy" operator insertions Hj (which are allowed to be 
built of both heavy and light fields) be divided into several non-intersecting groups (nu- 
merated by A) in such a way that the heavy momentum conservation holds separately 
within each group: 

J2Qj = 0(k), for each A. (6.23) 

In terms of the independent momenta Q : 

J2 Qj = 0, for each A. (6.24) 

For each group it is natural to introduce a separate test function F\(Q\) (where Q\ 
denotes the set of all Qj for j e A), and the corresponding integration dQ' x over the 
manifold described by (6.24). Now eqs.(6.9) and (6.10) are replaced by: 

G(F x ,L)=<TS^[F x *n x e c ^ >„, (6.25) 

where 

Fx * H x = [ dQ' x F x (Qx) U Hj(Q 3 + <&)■ (6-26) 

jex 

To use graphical representation, one introduces an F-vertex of the same type as de- 
scribed after (6.6) for each of the test functions F\. 



6.6 Generalized operator expansions 

The As-operation on a diagram T contributing to (6.25) will have the form 



As ° r = E Il A as°M X 



r/nfc 



(6.27) 

M=oo,F=0 
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(cf. (5.12)), where the heavy knots h are 1PI with respect to light lines and must contain 
at least one F-vertex or heavy line. Setting M = oo and F = is an expression of 
the fact that all the "heavy" elements (F-vertices and heavy lines) must be contained 
within the fa's. 

Applying As to (6.25), one obtains the following expression instead of (6.21): 

AsoG(F x ,£(<p,$)) 



T]jF x *Hxe 

AeA 



£( V ,*) 



xexpA as o[Te £ ^-l]}> , (6.28) 



=e<t{ n A a 

{A} Ae{A} 

where {A} denotes a splitting of the set of all A into non-intersecting subsets numerated 
by A, so that 

n n F A*ftA=n F A*ftA, (6.29) 

Ae{A}AeA A 

and the first summation on the r.h.s. of (6.28) runs over all such splittings {A}. 

The action of A as can be most easily described using the "fool-proof recipe pre- 
sented in subsect. 5.2: A^ nullifies unconnected diagrams as well as all those connected 
ones which result in meaningless expressions like 1/0 when formally Taylor-expanded 
in light parameters. The remaining diagrams are exactly 1PI with respect to light lines 
and get formally Taylor-expanded in the light parameters. 



Each A as results in a local functional (cf. (6.12)), so that (6.28) can be represented 



as 



AsoG({F A },£(^,$)) 



=e<t{ n 

{A} 



CA,a A ,barc({FA}AGA, M, C, {fc}jeA, e)ja A £ Si 

Ae{A} |_ «A v'eA 

X exp£ effj barc(^)} >o, (6.30) 

where each C is a coefficient function which is independent of fields and is at the same 
time a linear functional with respect to each F\. The dependence on qj is analytical. 

To obtain an expansion for the renormalized Green functions in an explicitly con- 
vergent form one should — as has been done in all the special cases considered above — 
"sandwich" the As-operation between R and R _1 (cf. (5.27)). Further reasoning is 
similar to that of subsects.4.7 and 5.2-3: one should first get rid of the i?-operation us- 
ing the generalized Zimmermann identities of subsect. 4.1, then apply the As-operation, 
and, finally, use the identities of subsect. 4.2 to extract the i?-operation and thus trans- 
form the expansion to an explicitly convergent form. In this way one arrives at the 
following final result: 

AsoRoG({F A },L(v?, $)+y?J) 



£Ro< T { n 

{A} Ae{A} 



E c A,a A ({FxjxeA, M, g, {^} jeA ) ja A [ E QJ 



x exp [L cS (ip) + ipj]} > , 



(6.31) 
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which should be compared with (6.17); L eS is given by (5.21) and has the form of 
(5.22). 

6.7 Example 

To clarify the above general expansion, consider the case of two groups of two currents 
each (the corresponding phenomenological problem is the deep inelastic scattering of 
two deeply virtual photons [13]). There are two heavy momenta Q\ and Q 2 , and two 
test functions, -Fab(Oi) an d Fcd(Q2)- The expression to be expanded is: 

J J dQ 1 dQ 2 F AB (Q 1 )F CB (Q 2 ) 

x Ro < T {Ja(Qi)Jb(-Qi + <h)MQ2)M-Q2 + g 2 )e^)+^} > , (6.32) 

where q are introduced to get rid of the non-connected components. Its expansion 
is as follows. There are two sets {A} over which the summation in (6.31) runs: 
{(A, B), (C,D)} and {A, B, C, D}. The first of them corresponds to that term in the 
expansion where the two pairs of currents are "shrunk" into (linear combinations of) 
local operators separately; the second corresponds to a single sum of local operators. 
The resulting expansion reads (cf. Fig. 6): 

J J dQ 1 dQ 2 F AB (Q 1 )F CD (Q 2 ) 
x EE c ab,o(Qi, M, g, fi)c C D,b(Q2, M, g, //) 

\ a b 

xR° < T {j a {qi)jb{q2) exp[L cff (v9) + y?J]} > 
+ E c ABCD,a(Qi,Q2,M,5(,/i)R° < T {j a (?i + ^) exp[L cff (ip) +ipJ}} >oj- (6.33) 

6.8 Contact terms in operator expansions 

The expansions that we have derived are valid for arbitrary test functions F(Q). For 
example, if there are two independent heavy momenta, Qi and Q 2 , then the expansion 
will in general contain terms proportional to S(Q 1 ), 5(Q 2 ), S(Q 1 + Q 2 ), S(Qi)S(Q 2 ) etc. 
(for a discussion of the role of such terms see [1]). Let us explain how such terms are 
generated by the As-operation. 

Within the context of subsect. 6.2 (no heavy masses etc.), consider the case corre- 
sponding to the product of two currents in (6.8): 



F*H = J dQ F(Q)H(Q + q)H(-Q). 



(6.34) 
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The non-trivial terms are generated from the following expression (cf. (6.12)): 



TF * He 1 



= Y, c n7( F , L )Uq), (6-35) 



where L is the interaction Lagrangian. 

There are two classes of terms generated by TF*He L that give non-zero result after 
application of A^. The first class contains the terms whose graphs remain connected 
after the F-vertex is deleted. Such terms can be represented as: 

^ J d D QF(Q)f(Q, Pl , . . . ,p n ,m,/i) x S(q + U(<p( Pj ) d D p 3 ), (6.36) 

3 

where the Feynman amplitude / depends on the masses m of the theory, contains one 
factor fi 2e per each loop, and the 5-function that expresses the momentum conservation 
is shown explicitly. The action of A as on (6.36) consists in Taylor-expanding / in masses 
and pj. The result has the form: 



A as o(6.36)= / d D QF(Q)J2fa(Q,ri 

J n 



x 



H 2e J V a (m,p 1 ,...,p n )5(q + J2Pj)I[( l P(Pj) dD Pj) 



(6.37) 



The square-bracketed term in (6.37) corresponds to j in (6.12), and the dependence of 
fa(Q,fi) on Q an d H is as follows: 



2 

m 



fa(Q,fi)=fa(Q,{-Q 2 -) )• (6.38) 

In the /-loop approximation, f a is a polynomial of order / in its second parameter, while 
its dependence on the first parameter is determined by power counting and covariance 
properties. 

The second class contains the heavy knots that consist of two parts connected only 
via the F- vertex (cf. Fig. 7). If there are no heavy masses, then there is only one 
such subgraph that is not nullified by the operation r — it consists of the two vertices 
corresponding to the "heavy" operators H (see (6.34)) connected to the F-vertex. The 
corresponding contribution to (6.35) has the following form (for definiteness we assume 
that H(x) = (p(x) 2 ): 



A as o[ J dQ F(Q) S(q + Q+p 1+ p 2 )5(-Q + p[ + p' 2 ) 



X(fi(p 1 )(p(p2)(p(p' 1 )'fi(p' 2 ) dp[ . . . dp' 2 \ 
r[dQF(Q)5(-Q+p' 1 +p' 2 )] 

X S(q +pi +p 2 +p[ + P2)v(Pi)v(P2)<p(Pi)<p(P2) dpi... dp 2 , (6.39) 
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where the two ^-functions in the first line are inherited from the two operators H (cf. 
(3.1)), and we have extracted a 5-function in the last line that depends only on the 
light momenta. We have deliberately not performed the integration over Q on the r.h.s. 
because r should Taylor-expand the square-bracketed expression in p\ and this can be 
done by applying r directly to the 5-function. The result has the form 

J2 I dQ[F(Q)5^(Q)]xj a (q), (6.40) 

a J 

and one sees that the "bare" coefficient functions in (6.12) receive 5-functional contribu- 
tions. Such contributions — "dressed" by the appropriate divergent renormalization fac- 
tors (cf. (6.18)) — are finally inherited by the coefficient functions of the short-distance 
expansion (6.17). It may be not quite obvious that the coefficient functions thus con- 
structed will be integrable around Q ~ 0. However, cancellations between different 
terms can be traced if one performs the expansion procedure explicitly starting from 
the As-operation for products of singular functions prior to performing the integra- 
tions of ^-functions described in subsect. 5.3. Alternatively, an interested reader can 
verify validity of our recipes by straightforward calculations in a simplest situation, 
e.g. within the model ip 3 in two dimensions for the currents H(x) = <p 2 (x) in one loop 
approximation. 



6.9 General linear restrictions on heavy momenta and paralo- 
cal operators 



Consider the following operator product instead of (6.10): 



F *H 



J dQ'F(Q) 



UH S (Q S + qj) 



(6.41) 



where the ^-functions are introduced to impose linear restrictions on the heavy mo- 
menta. With a special choice of c X j one reproduces the restrictions (6.24). Examples 
of more general restrictions are presented in Fig. 8. We are not aware of any phe- 
nomenological applications where such restrictions might emerge naturally. Therefore, 
we offer only a few comments concerning the most general case (6.41). 

First, it should be stressed that the expansion procedure as described by (6.11) 
remains fully correct and well-defined provided the "fool-proof recipe for A as from 
subsect. 5.6 is used. No simple graphical description for heavy knots exists, though. 

Second, the representation (6.12) should be replaced by a more general one: the 
sum on the r.h.s. may now contain operators that can not be interpreted as local 
operators and differ from (3.1) by additional 5-functions besides the one expressing 
momentum conservation. Such operator monomials can be called paralocal operators 
[3]. A simple example of a paralocal operator in position-space representation may be 
as follows: 
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but more complicated patterns are possible: note in this respect that a natural position- 
space representation for the multilocal operators which are a special case of paralocal 
ones, is 

Ji(zi)..J n (z n ), (6.43) 

where all J are local operators. 

The renormalization properties of paralocal operators can be studied most easily 
using the techniques of [15] (also [16]) and are similar to those of (multi-) local ones: a 
paralocal operator renormalizes via paralocal ones. Zimmermann identities (including 
the inverted ones) can also be generalized to this case. Therefore, expansions in explic- 
itly convergent form for the case of general non-natural restrictions on heavy momenta 
can be obtained in a manner completely similar to the case of natural restrictions. 



Conclusions 

We have developed a simple combinatorial technique for studying global exponen- 
tiation properties of the As-operation and obtained Euclidean asymptotic expansions 
for MS-renormalized Green functions of arbitrary local operators in arbitrary models. 
The expansions are true infinite asymptotic series that run in powers and logs of the 
expansion parameter. The obtained expansions exhibit perfect factorization of heavy 
and light dimensional parameters (masses and external momenta) which means e.g. 
that the coefficient functions of operator expansions are analytical in light masses and 
momenta. It should be stressed that the expansions are valid in models with massless 
particles like QCD, and are most convenient for practical calculations. Uniqueness of 
such expansions (cf. [1]) greatly facilitates study of their properties in gauge models. 
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Figure captions 

Fig. 1. A diagram (a) and its heavy knots (a)-(d). All the external momenta are 
light. The fat lines correspond to propagators of heavy particles. Various IR-subgraphs 
are shown with dashed lines in (b)-(d). 

Fig. 2. (a) An example of a diagram with a "light" tadpole subgraph (the notations 
are the same as in Fig. 1). The tadpole corresponds to a non-expandable global factor. 
The As-operation should not affect such factors, (b) and (c) are two IR-subgraphs. 
The contribution corresponding to the IR-subgraph shown in (b) is automatically set to 
zero by the operation r due to the properties of dimensional regularization, so that the 
tadpole can be factored out in the result produced by the As-operation, as expected. 

Fig. 3. A graphical representation of the test function corresponding to the heavy 
external momenta. The lines connecting the F- vertex with the rest of the diagram are 
heavy by definition. 

Fig. 4. If the wavy lines correspond to massless scalar particles then setting q = 
would result in an IR divergence in D = 4. 

Fig. 5. An example of kinematics of heavy momenta with linear restrictions, (a) 
differs from (b) by the restriction Q 3 = —Q2- (c) shows the corresponding position 
space picture: £, ( — > 0. 

Fig. 6. A diagrammatic illustration of (6.33). The two terms on the r.h.s. corre- 
spond to two ways of cutting off the heavy knots corresponding to the heavy external 
momenta. 

Fig. 7. Graphical representation of the heavy knots contributing to contact terms. 
Each blob corresponds to a subgraph that is 1PI with respect to light lines. If there 
are no heavy masses then each blob is just 1PI and is nullified by r unless it consists 
of a single vertex. 

Fig. 8. Examples of "non-natural" restrictions on heavy momenta, (a) differs from 
Fig. 6 by an additional restriction Q\ —Qi- 
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